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Two unsteady flows dominated by the occurrence of separation are simulated 
through the use of the discrete v’ortex model. The first of these is a sinusoidally 
oscillating flow about a circular cylinder at a Keulegan-Carpenter number of K= 10. 
The vorte.x model has been combined with the boundary' layer calculations and the 
positions of the separation and stagnation points, the evolution of the wake, the 
velocity and pressure distributions, and the instantaneous forces have been calculated 
and compared, whenever possible, with those obtained experimentally. The model has 
successfully simulated the occurrence of the transverse half Karman vortex street. The 
calculated positions of the vortices were found to be in good agreement with those 
obtained experimentally. The measured and calculated in-line forces and the 
dilTerential pressure distributions showed reasonably good agreement. 

The second simulation dealt with a rapidly decelerating How about a two- 
dimensional sharp-edged camber. An extensive study of the velocity field in the 
vicinity of the singular points led to the development of a novel method for the 
introduction of vorticity at variable time intervals. The measured and calculated 
characteristics of the flow, such as the evolution of the wake and the forces acting on 
the camber, were found to be in excellent agreement. Furthermore, the simulation 
provided a plausible explanation for the cause of parachute collapse, a phenomenon 
which gave impetus to the numerical and physical e.xperiments described herein. 
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I. DISCRETE VORTEX ANALYSIS OF SINUSOIDALLY OSCILLATING 

FLOW 

ABOUT CIRCULAR CYLINDERS 



A. INTRODUCTION 

The separated steady and unsteady flows about bluff bodies have been almost 
completely unyielding to both analysis and numerical simulation for a number of 
mathematical reasons and fundamental fluid dynamic phenomena. Separation gives 
rise to the formation of free shear layers which roll up into vortex rings or counter- 
rotating vortices. They, in turn, interact with each other, with the counter-sign vorticity 
generated at the base of the body, and with the motion of often unknown separation 
points. The wake becomes unsteady even for a steady ambient flow and the problem 
of the determination of the characteristics of the wake becomes coupled to the 
conditions prevailing upstream of the separation points. Evidently, viscosity modifies 
radically the inviscid flow, which, in this case, cannot serve even as a first 
approximation to the actual flow. The boundar}' layer equations are not applicable 
beyond the separation points and are, therefore, of limited use in bluff-body flow 
problems. 

Page and Johansen's pioneering experimental work (192S). Gerrard's (1966) 
vortex formation model, and Roshko's (1954) numerous contributions, followed by a 
large number of important papers, have provided extremely useful insights into the 
mechanism of vortex shedding. It became clear that a two-dimensional body immersed 
in a two-dimensional steady flow does not give rise to a two-dimensional steady wake 
and only a fraction (about 60% for a circular cylinder) of the original-circulation 
survives the vortex formation. It also became clear that bluff-body flows exhibiting 
separation, turbulence, and time-dependence are almost completely unyielding to both 
analysis and simulation even if the ambient flow is assumed to be time invariant. 

Many flows of practical interest are unsteady, i.e., the characteristics of the 
amxbient flow are time-dependent. In the past twenty years or so a large number of 
theoretical and experimental studies have been carried out. These dealt primarily with 
unseparated laminar flows, the early stages of impulsively started flow over plates and 
cylinders (numerical and experimental studies), and oscillating flows with zero or non- 
zero mean flow (on an infinite flat plate and over a cylinder with streaming flow, all 
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under laminar flow conditions) for the purpose of studying the eflects of flow 
unsteadiness on the transition mechanism and turbulence development (see e.g. 
Bradbur\’ et al. 1982). Ver\’ little has been attempted either theoretically or 
experimentally to analyze the wake-bouhdar\'-layer interaction in time-dependent flows 
(i.e., with unsteady ambient flow). 

The subject of separated time-dependent flow at large Reynolds numbers is lesser 
developed but of greater practical importance (particularly to marine related topics) 
relative to other classical component disciplines of fluid mechanics. 

A number of unsteady flow machines and their use in the investigation of 
unsteady turbulent boundarx’ layers have been reviewed and documented by Carr 
(1981). These included flat plate, diffuser, pipe, airfoil, and cascade flows. The results 
have shown that (i) the time-averaged mean velocity profile is almost always the same 
as the velocity profile that would occur in a steady flow having an equivalent mean 
e.xternal floxv velocity; (ii) the turbulent structure in the oscillating flow is not changed 
from the equivalent steady-state counter part; and (iii) the unsteady eflects are often 
confined to a thin layer near the wall, while the outer region of the boundarx' layer is 
not strongly affected. These conclusions, apparently valid for unsteady turbulent 
boundarx' layer floxvs. are not applicable to unsteady separated, turbulent, blulT-body 
floxvs. 

The separated unsteady floxv situations involving xx’ake return, as in the case of a 
sinusoidally oscillating floxv about a cylinder, or xx'ake retardation, as in the case of a 
decelerating parachute, are an order of magnitude more complex. 

In steady flow the position of the separation points is nearly stationarx’, except 
for small excursions of about 3 degrees (on a circular cylinder). Furthermore, the 
interference betxveen the vortices and the body is confined mostly to the vortex 
formation region. 

For oscillating floxx's the net effect of the shed vortices is txvofold. Firstly, their 
return to the body dramatically affects the boundarx’ layer, outer flow, pressure 
distribution, and the generation and survix’al rate of the new vorticity. Secondly, they 
not only give rise to additional separation points (during the early stages of the floxv 
reversal) but also strongly aflect the motion of the primarx- separation points. These 
effects are further compounded by the diffusion and decay of vortices and by the three- 
dimensional nature of the floxv (all of xvhich gix-e rise to cycle-to-cycle x’ariations, 
numerous floxv modes, etc.). The stronger and better correlated the returning vortices. 
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the sharper and more pronounced the changes are in the pressure distribution on the 
body and in the integrated quantities such as lift, drag, and inertia coefficients. 
Nevertheless, the increased correlation does not entirely eliminate the consequences of 
the stochastic variations in the motion of vortices. 

In periodic flow, the mobile separation points (when they are not fixed by sharp 
edges), undergo large excursions (as much as 120 degrees during a given cycle of 
oscillating Oow over a circular cylinder). This experimental fact renders the treatment 
of boundaiy layers on blulT-bodies subjected to periodic wake return extremely difl'icult. 
particularly when the state of the boundary layer changes during a given cycle. 
Furthermore, the classical criterion of separation for steady flow, i.e., the vanishing of 
skin friction on the body, is no longer valid for unsteady How. According to the MRS 
criterion (Moore 195S. Rott 1956. and Sears 1972), it is the simultaneous vanishing of 
the shear and velocity at a point within the boundary' layer that determines the 
separation point. Furthermore, the time rate of change of circulation is no longer 
siven bv dF dt = 0.5L' " as in steadv flow, where U is the outer flow velocity at 

^ . S " S 

separation, but by (0.5U^ — where c9^ is the speed of the separation point. It is 

clear from the foregoing that there is little hope of devising a satisfactorx' theoretical 
model before something is understood of the unsteady processes associated with the 
formation and reversal of the wake, spanwise coherence, the sensitive dependence of 
the motion of vortices on small changes in the previous conditions and on the nature 
of transition in oscillating flow about smooth and rough cylinders. 

B. NUMERICAL SOLUTION OF THE N.AVIER-STOKES EQUATIONS 

For steady ambient flow about blufl' bodies (mostly cylindrical and airfoils), the 
numerical studies based on the use of the steady or unsteady form of the Navier-Stokes 
equations and some suitable spatial and temporal differencing schemes are limited, out 
of necessity, to relatively low Reynolds numbers (less than about 1,000) (see e.g., 
Lecointe & Piquet (1984) for a finite difference solution and Gresho et al. (1984) for a 
finite element solution). The major obstacles to the application of either the finite 
difference or the finite element methods to higher Reynolds number laminar flow are 
stability, computation time, treatment of the boundaiy conditions, and accuracy. Even 
though many differencing schemes have been developed to overcome the instability 
problem (Roache 1976). maintaining stability continues to be a problem with 
increasing Reynolds number. The truncation errors decrease the apparent Reynolds 
number by introducing an unknown artificial viscosity. Even if the problems 
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associated with stability and truncation errors were to be resolved, the attempt to 
obtain solutions of the Navier-Stokes equations at higher Reynolds numbers are 
limited by a fundamental fluid dynamic phenomenon: the stability of the flow itself 
When the flow becomes turbulent either in the wake and or in the boundary layers, 
one needs a closure model for turbulence to solve the Reynolds equations for a time- 
dependent, three-dimensional, separated, turbulent flow' (even if the ambient flow is 
smooth and the bluff body is two dimensional). Clearly, the roots of the most serious 
problem in the solution of the Xavier-Stokes equations are buried in the physics of 
turbulence. The stability and truncation-error problems associated with the 
differencing schemes may be resolved but the problem of turbulence appears to 
transcend all efforts. 

The numerical solution of unsteady incompressible Xavier-Stokes equations in 
their vorticity-stream-function formulation has been investigated by numerous 
researchers through the use of various finite-difference techniques. These studies 
concern mostly the separated flow about circular cylinder and prisms at relatively low 
Reynolds numbers (see e.g.. Davis & Moore 1982). 

It appears that the e.xisting numerical methods cannot yet treat the high 
Reynolds number flows with sufficient accuracy for a number of reasons. The finite 
difference schemes require a very fine grid, a turbulent model, and a ver\’ large 
computer memoiy. It seems that the modelling of the turbulent stresses in the wake, 
particularly in time-dependent flows, will be the major source of difficulty in all future 
calculations. Whether or not it will ever be practical to apply the finite difference and 
finite element methods to high Reynolds number flows is unknown. The inherent 
difficulties are certainly significant enough to warrant exploring other solution 
methods. 

C. VORTEX METHODS 

1. Discrete Vortex Analysis 

Certain separated time-dependent flows may be simulated through the use of 
the discrete vortex model (DVM) (see e..g. Chorin 1973; Sarpkaya 1975). The free 
shear layers which emanate from the sides of the body are represented by an assembly 
of discrete vortices. It has not yet been proven that a continuous vortex sheet may be 
so discretized. Thus, attention is given here to large scale vortex structures rather than 
to small-scale instabilities resulting from the vortex interactions. Furthermore, only the 
most important features of the method, as it is applied to sinusoidally oscillating flow 
about a smooth cylinder are described. 
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The previous attempts to apply the DVM to oscillating flow had either gross 
simplifications or met with various difficulties. Ward and Dalton (1969) considered 
only the symmetric flow situation with fixed separation points. Stansby (1977, 1979, 
1981) fixed the separation points at ± 90 degrees and used the velocity of the nascent 
vortex rather than the velocity at the separation point to calculate the vortex strengths. 
This resulted in significantly less vorticity input and prevented the returning vortices 
from interacting freely with the boundary layers and separation points. Subsequently. 
Stansby and Dixon (1983) used a Lagrangian vortex scheme and replaced the body 
surface by a polygon of line segments. The strengths of the segments of vortex sheet 
needed to establish the zero-velocity condition on the surface were determined from the 
inverse of an influence matrix and the tangential velocities just inside the surface. 
Subsequently, the segments were replaced by one or more point vortices. Stansby's 
calculations for Keulegan-Carpenter number K = 10, (K = U^^T/D), failed to predict 
the transverse vortex street observed experimentally. 

Sawaragi and Nakamura (1979) determined the separation points using 
Schlichting's (1932) periodic boundary layer theory (valid only for K < < 1). They have 
not used the Kutta condition and incorrectly included an image vortex at the center of 
the cylinder. Finally, the calculations were performed for only three-quarters of a 
cycle, hardly enough time for the transient flow to develop into quasi-steady state. 
Kudo (1979, 1981) investigated the sinusoidally oscillating flow about a flat plate 
normal to the flow. The wake was assumed to remain symmetrical. Kudo's model 
used a Kutta condition, combined with a highly complicated force-and momentum-free 
nascent vortex placement scheme. Ikeda and Himeno (1981) studied the oscillating 
flow about a cylinder and a Lewis form. Separation points were assumed to be given 
by Schlichting's (1932) solution. As with Sawaragi and Nakamura (1979), they have 
incorrectlv retained the image vortices at the center of the evlinder. 

The discrete vortex model, as used in the present investigation, is relatively 
free from the arbitrary assumptions and inconsistencies noted above. However, the 
inclusion of the effect of turbulence in the boundary layers and in the wake remains 
unresolved. .As it will be noted shortly, the determination of the separation points 
requires the use of a separation criteria based on a laminar or turbulent flow 
separation. 
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2. Formulation of the Problem 



The complex velocity potential may be written as 



1 1 1 

W(z) = -U(t)(z + — ) + — ^ [ Ln (z-z^ ) - Ln(z - — ) j (1.1) 

2 n= 1 

where U(t) = sin (27tt T). and z^ are the strength and position of the n-th vortex. 
The velocities are normalized by and the distances by c. The complex velocity is 
given by dW;dz = — u + iv. 

The instantaneous force acting on the cylinder may be calculated either 
through the integration of pressure or through the use of the rate of change of impulse. 
It is relatively easy to show that the pressure on the cylinder is given by 



: m 



CL(t) \ " / Z Z \ 1 T 

P(z) = - 2p ^ COS0 - pM [— ^ r (( — 2 ) + — -2— — j- — pq-(z) ] 



( 1 .: 



z = re'® 



In terms of dimensional quantities, the use of the rate of change of impulse 



yields. 



D + iL = ip — y r (z 




+ 27tpC-L' 



(1.3) 



where D and L represent, respectively, the in-line and transverse forces. 

3. Specific Details of Creation and Convection of the Vortices 
The solution procedure employed was as follows: 

(1) The positions of the stagnations points at the upstream and downstream faces 
of the cylinder are calculated. For this purpose, the points at which the 
velocity is zero and changes its sign are located, starting from the most recent 
stagnation points. When the ambient flow reverses its direction (at the 
beginning of each cycle) the stagnation points switch their positions (the 
upstream one becomes downstream and vice versa). 
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(2) The positions of the primar}’ and secondar}' separation points are calculated 
through the use of one of the following methods: 

(a) Pohlhausen's Method: The velocity is calculated at one degree intervals 
along the cylinder, starting at the stagnation points. Then the positions 
of the separation points are determined through the use of Pohlhausen's 
technique (for details see e.g., Schlichting 1932). When the vortices 
returning to the cylinder cause irregularities in the velocity distribution 
along the cylinder (because of the disproportionately large influence of 
some of the point vortices), the Pohlhausen's method fails to predict a 
separation point. Under these circumstances, the maximum velocity 
criterion is used. 

(b) Maximum Velocity Criterion: The velocity distribution is calculated 
along the cylinder and the position of the maximum absolute velocity on 
each side of the cylinder is determined. Then the separation points are 
located at points where the tangential velocity is a certain time-averaged 
fraction of the maximum velocity. Further details of this fraction will be 
discussed later. 

(c) .Absolute Maximum Pressure Criterion: When the difference between 

the primarx’ separation angles calculated at times t^ and either 

through the use of Pohlhausen's technique or through the use of the 
ma.ximum velocity criterion, is larger than 8 degrees, then the positions 
of the maximum absolute pressure are used to calculate the separation 
points. For this purpose, all velocity peaks are calculated through the 
use of the tangential velocity distribution. Then the maximum velocity 
at which the maximum absolute pressure occurs is determined. Then the 
positions of the separation points are calculated using the maximum 
velocity criterion. In other words, the pressure distribution is used to 
locate the position of the maximum velocity nearest the true separation 
point. 

The foregoing, relatively time consuming, steps could have been eliminated 
through the use of a suitable numerical filter so as to remove the secondarx' 
oscillations imposed on the velocity distribution by a fexv point vortices in the 
vicinity of the cylinder. It was realized that the use of such a filter will not 
only consume more computer time but will also introduce uncertainties into 
the calculations. Thus, its use xx-as disregarded. 
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(3) The velocitv U is calculated at each separation point and the ratio U L' is 
determined. Then the cumulative average of this ratio is calculated for use in 
connection with the maximum velocity criterion to determine the position of 
the separation points whenever the Pohlhausen criterion fails. 

(4) The strength of the primary nascent vortex is calculated from. 

r„p= (»5|L',W-L-A)-At (1.4) 

where L'^cB^ accounts for the motion of the separation point relative to the 
flow velocity prevailing at the separation point. The strength of the secondary 
vortex is calculated from. 



r = 0.5 |U I U • At (1.5) 

os ' s‘ s ' ^ 

In other words, the effect of the relative motion of the secondary separation 
point on is ignored. 

(5) The nascent vortices are placed at a distance £ along the radial lines passing 
through the separation points, i.e., 3.1 = (1-^c) exp(i 6^). The value of c 

which satisfies the Kutta condition is given by (see e.g., Sarpkaya & Shoaff 
1979) 



1 + ( irj 2 HU.) ^ 

1- (irj/2KLy 



( 1 . 6 ) 



(6) The velocity at any point r<r* is calculated using the exact solution of the 
Navier-Stok.es equations for a single rectilinear viscous vortex (Lamb 1932). 
i.e., -u + iv = (r 2iTr) [1 — exp( — r",4vt)] e"’*^ and at any point r>r* using the 
point vortex relation (-u + iv = dw'dz). The radius r is measured from the 
center of the vortex and 0 is taken plus in the CCW direction from the x axis. 
The core radius r* at which the tangential velocity is maximum is given by 
r* = V f 26(vt) where t is the age of the vortex since its inception. Then the 
vortices are convected through the use of a simple Eulerian scheme for a time 
interval At = 0. 125c U T 

m 
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(7) The vortices which come closer than a distance of 0.05c to the cylinder surface 
as well as those which have a strength smaller than 0.005 are removed from 
the field. 

(S) Oppositely signed vortices are combined when their separation distance is less 
than 0.1. Furthermore, the secondar>’ vortices are combined at suitable 
intervals with the primar\' vortices to reduce the number of the vortices and 
hence the computation time. 

(9) The pressure and tangential velocity distributions and the force acting on the 
cylinder {both through the integration of pressure and the rate of change of 
impulse, given by Eq. (1.3)) are calculated. 

(10) Ai the end of the foregoing calculations, plots of the vortex positions, 
tangential velocity, total pressure, and stagnation and separation points are 
produced. 

Two additional details of the model need to be described; The introduction of 
asynmietiy and the circulation dissipation. Any flow started impulsively from rest 
remains symmetrical (at least in its numerical simulations) if not disturbed for a short 
time period at or near the beginning of its inception. Sarpkaya and Shoaff(1979) have 
investigated the methods of asymmetiy introduction and found that the displacement 
of the vortices in one shear layer for a short time interval is much more suitable than 
any of the methods previously used. In the present calculations, a similar method was 
used. The only diflerence was that asymmetiy was introduced during roughly the first 
quarter cycle of flow (0.1 < t T < 0.3) in order to allow the asymmetry' to take effect in a 
reasonable time period. 

The previous investigations by Sarpkaya & Shoafl'(1979) and others (e.g.. 
Kiya & Arie 1977) have shown that the incorporation of dissipation into a discrete 
vortex model resulted in a reduction of lift and drag force magnitudes, but did not 
significantly affect the flow kinematics. In the present work, sample calculations with 
or without dissipation have shown also that the results differ only in the magnitude of 
the force coefficients and that the kinematics of the flow is not affected by dissipation. 
It is in view of this fact that in the results presented herein it was decided to avoid a 
relatively subjective dissipation mechanism. 

4. Discussion of Results 

The numerical calculations were carried out for a Keulegan-Carpenter number 
of K = 10. This value of K was chosen primarily because of the fact that some of the 



most important phenomena take place in the range 8 < K < 15 (Sarpkaya 1985: Bearman 
1985; Williamson 1985). The most important of these phenomena is the occurrence of 
a transverse half Karman vortex street on one side of the cylinder. 

A series of flow visualization experiments at K = 10 was conducted in a water 
table. Figures 1.1 and 1.2 show a time sequence from a representative run. In these 
experiments the cylinder (D = 1.5 inches) was oscillated in the water table using a 
period of T = 3 seconds (Re = 12.800 and Re/K = P = 1280). The end of the 
plexiglass cylinder was within 1, 16 inch of the bottom of the water table, so that the 
end effects were minimal. It can be seen that the vortex shedding is indeed on the 
same side (i.e.. left) of the cylinder. The numerical simulation of this phenomenon 
through the use of the discrete vortex model became a challenge to numerous workers 
(Sawaragi & Nakamura 1979; Ikeda & Himeno 1981; Stansby 1977, 1979. 1981; and 
Stansby & Dixon 1982, 1983). However, none of these investigators has been able to 
simulate the half Karman vortex street for a number of reasons described in the 
Introduction. The present simulation has removed the shortcomings of the previous 
analyses, discovered physically and theoretically realistic methods to deal with the 
separation points, eliminated the ad-hoc assumptions, and produced the results 
presented herein. 

Figures 1.3 through 1.7 show, at times T* = t/T = 0.2, 0.4875, 0.6063, 
0.7625, and 0.8625 the position of the discrete vortices, the tangential velocity 
distribution and the total pressure distribution. The arrow in each circle shows the 
direction and the magnitude of the ambient velocity. The vortices of opposite 
circulation are shown with different symbols. These figures will be used to familiarize 
the reader with the essential elements of the analysis. Subsequently, the evolution of 
the transverse vortex street, characteristics of the shear layers, and the normalized drag 
and lift forces will be presented. 

.As the ambient velocity starts from rest and nears its maximum value (Fig. 
1.3), the vortices grow nearly symmetrically on the downstream side of the cylinder, 
reminiscent of an impulsively started flow about a cylinder. This is a consequence of 
the spiralling of the primary shear layers (SL-1 and SL-2), emanating from the primary 
separation points (SP-1 and SP-2). The primary vortices give rise to secondary’ 
boundary layers on the downstream side of the cylinder. These, in turn, separate and 
give rise to two additional shear layers (called the secondary shear layers, denoted by 
SL-3 and SL-4). As time increases and the ambient velocity decreases (Fig. 1.4), the 




Figure 1.1 Time sequence from flow visualization experiments (0.12 < T* < 1.00). 
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Figure 1.2 



Time sequence from flow visualization experiments (1.13 <T* < 1.75). 
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first vortex sheds from the upper right hand side of the cylinder. This leads to the 
rapid growth of the vortex in the lower right hand side. The tangential velocity around 
the cylinder is now primarily due to that induced by the vortices. The spikes in the 
velocity and pressure distributions arc a consequence of the disproportionately large 
influence of a few vortices in the vicinity of the cylinder. They are most prevalent at 
times of very' small ambient velocity. 

For T* > 0.5. the ambient flow reverses its direction, the flow separates on 
the left side of the cylinder, and the previously shed vortices are convected towards left 
(Figs. 1.5 and 1.6) partly due to the ambient flow and partly due to the mutual 
induction of all the vortices and their images. The reversal of the flow has two major 
consequences. First, the convection of the previously shed vortex over the shoulder of 
the cylinder precipitates earlier separation and establishes, by its sense of rotation, a 
preferred position for the next dominant vortex. Second, the reversal of flow gives rise 
to additional primary and secondary separation points (Fig. 1.6). However, at a given 
time there are at most four shear layers, i.e., some of the separation points disappear 
and some new ones come into existence as the flow reverses. Subsequently, the 
vortical structure at the upper left side of the cylinder moves further downstream and 
the vortex on the lower left side of the cylinder (Fig. 1.7) grows rapidly (in a manner 
similar to that shown in Fig. 1.4). It begins to move towards the top of the cylinder, 
partly due to the mutual induction velocity of the aforementioned vortical structure 
and partly due to that of its image. The events just described more or less repeat 
themselves in the subsequent cycles (see Figs. A. 1 through .A..17 in Appendix ,A.). More 
or less, because no two cycles can be expected to be exactly alike due to the ever 
increasing number of vortices or. in other words, due to the transient nature of the 
flow (for an experimental confirmation of this fact see Sarpkaya 19S6). Many more 
cycles of calculations will have to be carried out in order to reach a quasi-periodic 
state. However, this is nearly impossible and hardly necessary, partly because of the 
computer-time limitations and partly because the calculations over a period of three 
cycles are more than sufficient to delineate the fundamental characteristics of the flow. 
The reason for the latter is provided by the mechanism of vortex shedding during each 
half cycle of flow'. As noted in the course of the discussion of Figs. 1.3 through 1.7, 
the flow' in each half cycle behaves as if it started more or less anew, the single 
convected vortex over the shoulder of the cylinder precipitating earlier separation and 
dictating the side from which the next dominant vortex is to be shed. In numerical 
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Figure 1.3 Position of vortices, velocity and pressure distributions at T* = 0.200. 
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Figure 1.4 Position of vonices, velocity and pressure distributions at T"" = 0.4875. 
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Figure 1.5 Position of vortices, velocity and pressure distributions at T* = 0.6063. 
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Figure 1.6 Position of vortices, velocity and pressure distributions at T* = 0.7625. 
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Figure 1.7 Position of vortices, velocity and pressure distributions at T* = 0.8625. 
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calculations the vortices do not dissipate and, therefore, begin to cover the entire flow 
field. The real vortices decay rapidly under the inlluence of viscosity (no vortex in the 
far field survives more than a few cycles). The purpose of the numerical calculations 
was not to simulate the effects of viscosity but rather to show that the behavior of flow 
(e.g., the evolution of the transverse vortex street) is primarily an inviscid phenomenon, 
dependent on the relative amplitude of the flow or the Keulegan-Carpenter number. 
The results show that the effect of viscosity is important in bringing about the 
separation and in decaying the vortices but not in bringing about the particular vortical 
structure. 

The vortex motion about the cylinder is better illustrated through the use of 
the vector plots of the velocity field. Figures 1.8 through 1.20 show the velocity field 
at suitable time intervals over one and one-half cycles, starting at T* = 0.5. At the 
end of the first half cycle (Fig. 1.8), two large vortices are separated by a strong current 
between them. As the flow reverses (Fig. 1.9) and the ambient velocity reaches a value 
of IL'l = 0.707, the vortex adjacent to the cylinder splits into two unequal parts as it 
moves to the left, partly over the top and partly over the bottom of the cylinder. The 
vortex to the right of the current in Fig. 1.8 is simply convected towards the left of the 
cylinder. At T* = 0.75 (|L'l = 1), the vortices convected over the top of the cylinder 
give rise to earlier separation on the upper left side of the cylinder. Furthermore, the 
new vortex growing on the lower left side of the cylinder is pulled towards the top of 
the cylinder. Figures 1.11 (1U| = 0.707) and 1.12 (1U| = 0) show the development of 
two large vortices, again separated by a strong current (compare Figs. 1.8 and 1.12). 
Figure 1.13 (IL'l = 0.707) shows that the vortex adjacent to the left side of the cylinder 
remains nearly intact and is convected towards right over the top of the cylinder. 
Figures 1.14 (T- = 1.25, jL'l = 1) and 1.14 (T* = 1.375, |L'l = 0.707), show that the 
separation of flow on the right side of the cylinder gives rise to new vortices, in a 
manner similar to that shown in Fig. 1.10. Subsequently, this leads to two large 
vortices (Fig. 1.16), separated by a current as in Fig. 1.8. Evidently, Figs. 1.8 and 1.16 
are not expected to be identical in an ever evolving unsteady flow. Nevertheless, the 
essential features of the two flows (separated by a full cycle) are quite comparable. 
The secondary vortical structures seen in Fig. 1.16 are a consequence of the inviscid 
nature of the calculations. .As noted earlier, they do not survive the effects of viscosity 
and turbulence in physical experiments. Figures 1.17 through 1.20 show the evolution 
of the flow for an additional half cycle. Clearly, the flow pattern repeats itself with 
minor differences (compare Figs. 1.8 and 1.16 and Figs. 1.12 and 1.20). 



It is evident from the foregoing that the reason for the occurrence of the 
transverse vortex street in the range 8<K<I3 (corresponding to the amplitude-to- 
diameter ratios of A/D = 1.27 to 2) is, as expected, the formation of two asymmetric 
vortices during a given half cycle and’ the effect of the returning vortices on the 
formation of the new vortices. 

The vortical structures and the transverse current shown in Figs. 1.8 through 
1.20 cannot be directly compared at an arbitrary time T* with those shown in Fig. 1.1 
since the streaklines of an oscillating flow about a stationary cylinder are not the same 
as those of an oscillating cylinder in a fluid otherwise at rest. Nevertheless, a 
comparison of Figs. I.ID and I.IL with Figs. 1.12 and 1.20 (at times of zero ambient 
velocity or zero cylinder velocity) has shown that the radial and the angular positions 
of the vortices are predicted with surprising accuracy. In Figs. I.ID and I.IL the 
vortex to the left of the current is located at r/c = 5.45 and 0 (from the vertical) = 17 
degrees. The position of the same vortex in Figs. 1.12 and 1.20 is r/c = 5.42 and 0 
(from the x axis) = 15 degrees. It should be noted in passing that the side from which 
the dominant vortex sheds may become switched by the action of random disturbances 
in the flow or by stopping and restarting the flow. In the numerical calculations the 
direction of the transverse flow depends on the side at which the asymmetry is 
introduced into the flow. Had the asymmetry been introduced on the lower side of the 
cylinder, the transverse flow would have occurred on the lower side of the cylinder. 

The creation, shedding and backward convection of the dominant vortex in a 
flow' field, w'here there are only a few vortices, strongly affect the pressure distribution 
about the cylinder. Thus, it w'ould have been desirable to compare the measured and 
predicted pressure distributions. However, no such data is available at K = 10 or at 
any other K value. Wilson (1983) measured the differential pressure between two 
diametrically opposed points on a cylinder and provided differential pressure 
distributions at suitable times. Calculations have been performed with the numerical 
model to obtain similar differential pressure data. Figure 1.21 shows two comparisons 
of the measured and calculated differential pressure distributions. Even though the 
magnitudes differ somewhat, the general shape of the two Cp curves are remarkably 
similar. The reason for the difference in their magnitude is directly related to the 
strength of the vortices. Previous investigations (see e.g., Sarpkaya 1973) have shown 
that the calculated vortex strengths are somewhat larger than those obtained 
experimentally. An additional loss of vorticity could have been introduced into the 
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model in a manner similar to that used by Sarpkaya and ShoalT (1979) in order to bring 
the strengths of the shed vortices, and thereby the differential pressure distributions, 
into closer agreement with those obtained experimentally. However, it was decided to 
avoid a relatively subjective dissipation mechanism. It appears that the kinematics of 
the flow does not strongly depend on the strength of the shed vortices. This is in 
conformity with the previous applications of the discrete vortex model (Sarpkaya & 
Shoafi'1979). 

Figures 1.22 and 1.23 show the calculated drag and lift coefficients as a 
function of time. The effect of the passage of a vortex over the shoulder of the 
cylinder is exhibited in these figures in two ways. First, it brings sharp changes in Cd 
near its maximum. Second, it gives rise to a net transverse force on the cylinder. 
Evidently, the lift force is ver\' sensitive to the flow asymmetrx’ and depends strongly on 
the proximity to the cylinder of the vortex passing over the cylinder. The calculations 
will have to be carried out over many more cycles in order to obtain a root-mean- 
square value of the lift coefficient for comparison with that obtained experimentally. 
This has not been done in the present study for the reasons cited earlier. 

■A comparison of the measured and calculated drag force is shown in Fig. 1.24 
The agreement is not as good as expected primarily due to the fact that the strength of 
the calculated vortices is somewhat larger than those encountered in the experiments. 
A 5 to 10 percent decrease in circulation could have brought the measured and 
calculated forces into much closer agreement. As noted earlier, however, this was not 
the purpose of the present investigation. 

The attention will now be directed to the motion of the stagnation and 
separation points and to the vorticity fed into the shear layers. Figures 1.25 and 1.26 
show the motion of the upstream (primaiy) and downstream (secondan,') stagnation 
points, respectively. Here upstream means the side of the cylinder facing the ambient 
How. The primary and secondary stagnation points are approximately ISO degrees 
apart, as would be expected. For example, at T* = 0.5 the primarx’ stagnation point 
jumps to nearly zero and the secondarx' stagnation point jumps to approximately 180 
degrees. The spurious oscillations in Figs. 1.25 and 1.26, particularly in the position of 
the secondarx' stagnation point, are due to the disproportionately large influence of the 
proximity to the cylinder of a few point vortices. 

Figures 1.27 and 1.28 show the motion of the primarx' and secondary 
separation points (see also Fig. 1.3). .Aside from the spurious secondary oscillations 



34 



(due to the reasons cited earlier), Fig. 1.27 shows in general that the primary' separation 
angle increases rapidly as the velocity of the ambient flow increases from zero (e.g., at 
T* = 1. see Fig. 1.12) to an average value of about 110 degrees at T* = 1.5 (see Fig. 
1.16). Then it decreases rapidly to zero as T* approaches 2 (see Fig. 1.20). .A, 
comparison of Figs. 1.27 and 1.28 shows that there is a phase shift of one half of a 
cycle between the positions of the primary and secondary separation points, as would 
be e.xpected. There is, at present, no separation point data for comparison with those 
obtained experimentally. 

Finally, Fig. 1.29 shows the vorticity fed into each shear layer. The first and 
second shear layers have relatively large and oppositely-signed vorticity. They are not 
each other's mirror image due to the asymmetrx' of the flow. The vorticity in the third 
and fourth shear layers is about 25 to 30 percent of that in the primarx' shear layers. 
Clearly, the vorticity fed into the primarv' shear layers (particularly that fed into the 
first shear layer) has a sinusoidal character to it. This is because the vorticity is related 
to the ambient flow velocity (sinot) through the velocity at the separation point. 
The flow variables such as vorticity cannot be measured directly. Thus, the accuracy 
of its magnitude and time-dependence can only be inferred indirectly through the 
comparison of the measured and calculated forces and pressure distributions. It 
appears, on the basis of such comparisons that, the calculated vorticity is about 10 
percent larger than that prevailing in physical experiments. The reasons for this are 
partly the diffusion of vorticity by viscosity and turbulence and partly the three- 
dimensional nature of the physical experiments. 

D. CONCLUDING REMARKS ON OSCILLATING FLOW ABOUT 

CYLINDERS 

The discrete vortex model used in the present investigation avoids many of the 
problems associated with the Eulerian finite difference and finite element methods but 
is subject to some of its own. Chief among these is the excessive computer time needed 
for the calculation of the convection of vortices and the difficulty to account for the 
effects of viscosity and turbulence. 

In the present simulation the discrete vortex model has been combined with the 
boundary layer calculations and the positions of the separation and stagnation points 
have been calculated for a Keulegan-Carpenter number of K = 10. The results have 
accurately predicted the formation of a half Karman vortex street in the transverse 
direction. The calculated positions of the vortices were found to be in good agreement 



with those obtained experimentally. The measured and calculated drag force and the 
difTerential pressure distributions showed reasonably good agreement. The measured 
in-line force was somewhat smaller primarily due to the fact that the vorticity fed into 
the shear layers has not been artificially -reduced to bring the measured and calculated 
forces into closer agreement. The results have also shown that the effect of the 
backward convection of a large vortex over one side of the cylinder is indeed ver>’ 
pronounced on all the measured and computed characteristics of the flow. This is one 
of the most important reasons as to why the Morison's equation (see e.g., Sarpkaya & 
Isaacson I9S1) fails to represent the in-line force acting on the cylinder with reasonable 
accuracy, particularly in the range 8 < K < 13. 
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Figure 1.8 Flow field at T*= 0.500. 
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Figure 1.9 Flow field at T*= 0.625. 
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Figure 1.10 Flow field at T*= 0.750. 
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Figure 1.11 Flow field at T*= 0.875. 
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Figure 1.12 Flow field at T* = 1.00. 
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Figure 1.13 Flow field at T*= 1.125. 
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Figure 1.14 Flow field at T*= 1.250. 



43 



X/C 




Figure 1.15 Flow field at T*= 1.375. 
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Figure 1.16 Flow field at T* = 1.500. 
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Figure 1.17 Flow field at T*= 1.625. 
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Figure 1.18 Flow field at T*= 1.750. 
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Figure 1.19 Flow field at T* = 1.875. 
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Figure 1.20 Flow field at T*= 2.00. 
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Figure 1.21 Comparison of numerical and experimental differential pressures. 
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Figure 1.22 The calculated drag coefficient. 
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Figure 1.23 The calculated lift coefficient. 
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Figure 1.24 Measured and calculated drag force. 
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Figure 1.25 Upstream stagnation point. 
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Figure 1.26 Downstream stagnation point. 
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Figure 1.27 Primary separation points. 
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Figure 1.28 Secondary separation points. 
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Figure 1.29 Vorticity fed into the flow field. 
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II. DISCRETE VORTEX ANALYSIS OF UNSTEADY FLOW 
ABOUT CAMBERED PLATES 



A. INTRODUCTION 

The determination of the deployment sequence of an axisymmetric porous 
parachute and the unsteady aerodynamic loads acting on it present a very' complex 
coupled problem. The development of an analytical or numerical model which takes 
into account the effects of porosity, gaps, and variable opening schemes would allow 
numerical e.xpcriments on a large class of parachutes, reduce the number of expensive 
field tests to a few judiciously selected ones, and enable the designer to calculate the 
time history' of the fall of the parachute and the strength required to survive the 
aerodynamic loads. However, the development of such a model is hampered by a 
number of difficulties. 

The previous models for parachute loads are based by and large on empirical 
assumptions (see e.g.. Heinrich and Saari 1987; Mcwey 1972). They rely on the 
observation that families of parachutes open in a characteristic length and seem to 
have aerodynamic properties that relate well to the projected area of the parachute. 
The apparent mass is assumed to be a function of the projected area only and is not a 
function of the prevailing flow characteristics. The vortex sheet analysis w'as used by 
Klimas (1977) to derive the acceleration-independent apparent mass coefficient for 
arbitrary-shaped axisymmetric surfaces. Muramoto and Garrard (1984) used a 
continuous-source model to predict the steady-state drag of ribbon parachutes. The 
analyses did not, how'ever, deal with the evolution of the unsteady wake and its 
interaction with the canopy. 

It is in view' of the foregoing that a fundamental study of the separated time- 
dependent flow about two-dimensional rigid cambered plates was undertaken. Clearly, 
the flow' about a rigid cambered plate is considerably simpler than that about a porous, 
a.xisymmetric, and flexible parachute and the results, regardless of the degree of their 
agreement with corresponding experiments, may not have direct relevance to the 
practical problem under consideration. But the object of this investigation was the 
understanding of the evolution of the w'ake under controlled conditions rather than to 
provide a design tool. It is hoped that an investigation of this type will reveal the 
underlying physics of the phenomenon (particularly that of the parachute collapse). 
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help to interpret the full-scale results and will provide inspiration for the development 
of more general vortex models with which the dynamics of axisymmetric, porous, and 
flexible parachute canopies can be investigated. Efforts directed towards the 
development of a general numerical model, driven by the ever-present pressures of 
practical considerations, are deemed somewhat premature. Such efforts will have to 
face not only the problem itself but also the deficiencies of the vortex models and 
attempt to address to both of them simultaneously. The model presented herein 
removes the ambiguities associated with the use of the discrete vortex model and 
provides results which are in excellent agreement with those obtained e.xperimentally. 
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III. ANALYSIS 



A. TRANSFORMATIONS AND THE COMPLEX VELOCITY POTENTIAL 

The calculation of the velocity of any one of the vortices and the force acting on 
the body requires a conformal transformation (in which the camber becomes a circle), 
a complex-velocity potential representing the vortices, their images, and the two- 
dimensional irrotational flow around the body, and the use of the generalized Blasius 
theorem. 

The flow in the circle plane may be transformed to that about a cambered plate 
through the use of two successive transformations, one from plane to the plane 
and the other from the plan to the z plane. These are given by (see Fig. 3.1) 
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7 




2ib 




Figure 3.1 Circle and physical planes. 
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and ^ + m 



(3.1) 



z = i;" 







Combining the two, one has a direct transformation from the plane to the z 
plane as 



Z = £^+m — 

(, + m 



(3.2) 



It is easy to show that the camber in the z plane is a circular arc. 

The y-axis in the z plane passes through the tips of the camber. It is 
advantageous to locate the origin of the coordinate axes at the geometric center of the 
camber, i.e., at the center of the circle part of which represents the camber. This is 
easily accomplished by shifting the origin of the coordinate axes by 



where z^ is the x coordinate of the origin of the circle in the z plan. Thus, one has 
b" 

which transforms the circle in Fig. 3.1a to the physical plane in Fig. 3. Id. Table 1 
summarizes the relationship between m, z', the included angle of the camber, b, and 
the radius of the camber. 



TABLE 1 

SU.VI.VIARY OF THE PAR.A.METRIC RELATIONSHIP 



m 


✓ 


2a 


b 


R= l,m 


cos60 = .5 


-1 


120 


.866 




cos45 = .707 


0. 


ISO 


.707 


V2 


cos30 = .866 


1,V3 


240 


.5 


2V3 
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The complex potential function W in the circle plane (see Fig. 3.1a) which 
describes a uniform flow U (assumed to be time-dependent) with a doublet at the 
origin to simulate the cylinder, clockwise-rotating vortices (called q-vortices), Fj^p 
counter-clockwise rotating vortices (called p- vortices), and the images of all the p-and 
q-vortices in the circle plane may be written as 




in which Fj^p and i^^.p represent respectively the strength and location of the k-th p- 
vortex, Fj.^ and the strength and location of the k-th q-vortex, and c the radius of 
the cylinder; an overbar indicates a complex conjugate. The need for the separate 
identification of the p-and q-vortices and for the singling out of one of the vortices in 
each shear layer (namely F^p and F^^. nascent vortices) will become apparent later. 

B. COMPLEX VELOCITIES OF VORTICES 

The convection of the vortices and the calculation of the forces acting on the 
body require the evaluation of the velocities at the vortex centers. For the velocities in 
the circle plane this reduces to subtracting from Eq. (3.5) the complex potential 
corresponding to the vortex for which the velocity components are to be determined 
and evaluating the derivative of the remaining terms at To determine the 

velocities in the physical plane, however, one has to subtract (iFj^;27t) Ln(z — Zj^) from 
Eq. (3.5) or, in terms of , the terms (see e.g., Sarpkaya 1967, 1975) 



iF, 



IK 



Ln - C) + 



iF. 



In 



Ln 1 1 + 



(C+ni)(C +m) 



(3.6) 
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It should be noted that the first term in Eq. (3.6) is the complex function 
corresponding to the k-th vortex in the plane. The second term appears merely as a 
consequence of the transformation used. 

The above procedure may be generalized as follows. Consider the potential 
function for a single vortex in the physical plane and ignore, for the time being, the 
multiplier in front of the logarithmic term (i.e., irj^/27t). Then one has 

Ln(z - z^,) = Ln[ 1(0 ~ fig ] with z = flQ (3.7) 

equation (3.7) may be written as 

Ln (z - z^) = Ln (i; - g + Ln (3.8) 



Evidently, the first term on the right hand side of Eq. (3.8) represents the vortex in the 
circle plane. Let us now examine the derivative of the second term with respect to z. 
One has. 



dW 
— P 



d 



dz 




Ao-fig 



fio-gg y 



dz 



(3.9) 



where dz, di^ = f (Q. In the neighborhood of the function f (Q may be expanded as. 



fio - fig = (i;-gf''{g 



(i; - g* 



ag 



(3.10) 



Thus, one has 



d\V 

- P T 1 TVJ 


1 '2 f (g 


1 


' — LIIU 

dz 


1 

1 1 
1 





(3.11) 



or 
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(3.12) 



dz 2r- {Q 



Thus, the complex velocity in the physical plane reduces to 



ir 



- u + iv = — [ W(L) - Ln (i; - Q ] 

y 2:t f(0 



ir f(Q 



^ = C 



(3.13) 



in which for a p-vortex 






b2 



(^kp+ 



and 



2b- 



the last term in Eq. (3.13) reduces to 



(3.14) 



(3.15) 



_ (~b-)(C^^^ + m) 

2;t [(C^p + m)- 4- b^ ]2 



(3.16) 



This result could have been deduced directly from Eq. (3.6). However, the 
generalization of the method enables one to apply Eq. (3.13) to any vortex for any 
transformation between the circle and the physical plane. 



C. KUTTA CONDITION 

The fact that the flow separates tangentially with a finite velocity at the edges of 
the plate (Kutta condition) may be expressed by requiring 



dW 

— = 0 at C = - m ±ib (3.17) 
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Thus, inserting Eq. (3.5) in Eq. (3.17) one has 



2" \?r^op 




2” -Co, 




m 



+y- 









-U(l- ^) =0. 




Equation (3.18) may be decomposed into two parts as 









2 ^ 




+ (-^0 + i ^ o ) = 0 - 



(3.18) 



(3.19) 



where the terms containing the strength of the nascent vortices represent the velocity 
induced at the tip of the camber by the nascent vortices and the term in parenthesis the 
velocity at the tip due to all other vortices (and their images), the doublet at the center 
of the circle in the plane and the ambient velocity. 

Equation (3.19) represents two coupled equations for the strengths and positions 
of the nascent vortices. Thus, the solution of the said quantities does, in general, 
require an iteration. However, this iteration may be avoided by noting that the 
velocity induced by a nascent vortex at the opposite tip is very small and certainly 
negligible. Thus, Eq. (3.19) for one of the nascent vortices may be reduced to 




2:t 








(3.20) 
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A similar expression may be WTitten for the other nascent vortex. The use of the Kutta 
condition, as expressed by Eq. (3.20), will be further explained following the discussion 
of the tip velocity. It suffices to note that all nascent vortices satisfying the Kutta 
condition do not yield either the same tip velocity or the same velocity distribution in 
the neighborhood of the tip. There are, in fact, certain preferred positions for the 
nascent vortices which yield physically realistic velocity distributions near the tips of 
the cambered plate. These nascent vortex positions will be discussed later. 

D. TIP VELOCITY 

According to the Kutta condition the tangential velocity at the tip is finite. The 
purpose of the following is to determine this finite velocity. It may be determined 
either through the use of fHopital's rule or through the use of a more general 
expression which is valid for all other transformations. 

The velocity at the tip is given by 



dW_ dW di; 
dz di^ dz 



at z = z ± 2ib 

t O 



(3.21) 



For an arbitrary point z, Eq. (3.4) yields. 



dC 1 ^/TF 

dz 2 2Vz “ 



(3.22) 



In general, one may write Eq. (3.22) as 



di; 1 



(z - z^) 



dz 2 2v(z — z^)" + 4b* 



(3.23) 



or multiplying both sides with J{z — z^^), one has 



fiz) = (z-z )^^2 ^ 
dz 



1 (z - z^) 

— (z -z,)‘^“ ± — 

7 F _ 



Vz -(Z^± 2ib) 



(3.24) 
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Thus, as ± 2ib, one has 



f(z,) = (z 




dz 




(3.25) 



or 



£l = ^ 2r-(z,) 

dz ~ 2 v'T^) " C - 

Expanding dW/dC in the neighborhood of one has 
dW // 

— = (^ - g w (g + .... 

QL, ‘ ^ 



(3.26) 



(3.27) 



Combining Eqs. (3.26) and (3.27) one finally has 



dWi d-W , 



dz I du' ‘ 



(3.2S) 



2 = Z. 



Noting that for the case under consideration f" (zj = ib/4, one has 



dW 

dz 



z = z. 



d‘W ib 



(3.29) 



Equation (3.29) yields the desired finite tip velocities. It is easy to show that it may be 
obtained directly from Eq. (3.21) through the use of I'l-Iopital's rule. 



E. TIME DEPENDENT-FORCES 

The force acting on the body in the physical plane may be calculated either 
through the use of the pressure distribution or through the use of the rate of change of 
impulse. 
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Bernoulli's equation for unsteady flow is given by 




V 



1 \ - 



) - ( 




+ 




I av 

- j — ds = f(t) 
1 



(3.30) 



where the indices indicate two points on the body in the physical plane. Since there is 
no pressure drop across the shear layer and since the integral term in Eq. (3.30) is zero 
at the tip (i.e., ds = 0), one has 



f(t) = 





(3.31) 



where V^j and V ^2 represent the tangental velocities on the upstream and downstream 
faces of the tip. It is important to note that f (t) in Eq. (3.31) is also the time rate of 
change of circulation, i.e., the rate at which vorticity is shed into the wake from the tip 
of the cambered plate. 

The normalized form of Bernoulli's equation between any two points m and n 
then becomes 







+ 






+ 




V 




(3.32) 



The integration of the differential pressure between the upstream and downstream faces 
of the camber yields the force components in the x and y directions, i.e., the drag and 
lift forces. 

The force acting on the body can also be calculated through the rate of change of 
impulse. It is given by 

F - 4Hpc2u (1 - r„ (z„ - z„,) I (3.33) 
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which mav be written as 



+ i ^L" 



F cU c m* 

^ = 2n(— ^)(— )(1 - -^) 



2pLVb 



U 



b- 



2c‘ 






r_JL [f(^) 



2b 5{U t;c) c c 






(3.34) 



in which is the reference velocity; U, the rate of deceleration of flow and z=f(Cj^), 
i.e., the transformation given by Eq. (3.4). Equation (3.34) may also be deduced 
directly from the generalized Blasius equation. It is important to note that the force 
calculated from Eq. (3.34) includes the effect of the rate of change of circulation 
between two successive time steps. Thus, it may be smaller or larger (depending on the 
sign of T) than the force calculated through the integration of the instantaneous 
differential pressure Eq. (3.32). This is because of the fact that the instantaneous 
pressure depends only on the prevailing flow conditions and does not account for the 
rate of change of total circulation between successive time steps. In the calculations to 
follow L’^ and c are taken as unity for sake of simplicity. 

F. METHOD OF CALCULATION 

The methods used in the past in the determination of the vorticity flux from 
sharp-edged bodies may be roughly classified into two broad categories. The first of 
these involves the use of variable nascent vortex positions (see e.g., Sarpkaya 1968, 
1975) and the second, the use of fixed nascent vortex positions (see e.g., Clements 
1973-1975). 

The method of fixed positions involves the selection of a suitable fixed point in 
the flow near the separation point and the use of the velocity at that point to 
calculate the rate at which vorticity is shed into the wake from 



5F 

5t 




(3.35) 



In this method the positions of the nascent vortices are the crucial parameters. The 
previous applications of this method did not examine the effect of the position of the 
nascent vortices on the velocity distribution in the neighborhood of the separation 
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point. Only the distance of the fixed point to the body was varied and bracketed 

between two subjective limits by comparing the calculated results with those obtained 

experimentally. In this method no interaction is allowed between the shed vortices and 
the amplitude of oscillation of the point or the time of appearance of the nascent 
vortices. Furthermore, the time interval is chosen more or less arbitrarily (Kiya and 
Arie 1977) (repeating a few calculations with a single program with only the time step 
changed and also by referring to the results of the previous investigations). Thus, the 

velocities at the outer edges of the shear layers are only indirectly related to the 

strength of the nascent vortices and the fixed time interval. Evidently, the velocities in 
the inner and outer edges of the shear layers, the time interval, the strength and 
position of the nascent vortices, and the Kutta condition are interdependent and that 
both the position of the nascent vortices and the time interval cannot be chosen 
arbitrarily, even if they are chosen judiciously on the basis of previous experience and 
trial calculations. 

Sarpkaya (1975) used the method of variable nascent vortex positions and 
determined the rate of shedding of vorticity from the relation 



dr _ \ , , 

^ “ T 



(3.36) 



where is interpreted as the velocity in the shear layers calculated by using the 
average of the transport velocities of the first four vortices in each shear layer. The 
positions of the nascent vortices are chosen so as to satisfy the Kutta condition at the 
edges of the body and thus they can move slightly with time. Thus, this method 
simulates in a satisfactory manner the mechanism of feedback from wake fluctuations 
to the fluctuations in the rate of circulation. The number of disposable parameters is 
reduced to a minimum and in this sense this method is superior to the method of fixed 
positions. However, the use of the average of the transport velocities of the first four 
vortices remained questionable. 

It was often assumed that the vorticity flux could not be calculated, at each time 
interval, as it is applied to sharp-edged bodies, through the use of the mathematically 
finite velocity occurring at the sharp edges of the body. This assumption was based on 
the fact that the separation points are singularities of the transformation used and the 
numerical procedures may not be stable. 
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It is on the basis of the foregoing that an original study was undertaken to 
establish once and for all a method whereby the nascent vortices may be introduced 
into the flow without any ambiguities. The method fmally arrived at will be explained 
through the use of a series of figures and velocity plots. 

Figures 3.2a and 3.2b show the tip region in the circle and physical planes, 
respectively. The regions A and B in Fig. 3.2a were discretized through the use of a 
suitable grid and a single vortex was placed at a grid point. The strength of the vortex 
was determined from the Kutta condition Eq. (3.20). Then the velocity normal to the 
radial line OZ in the physical plane (Fig. 3.2b) was calculated in the vicinity of the tip 
through the use of the complex velocity potential. 




Figure 3.2 Tip region in the circle and physical planes. 

Placing the vortex along the radial line OM (in the plane) yields a single valued 
tip velocity independent of the strength and the position of the vortex and dependent 
only on the plate geometry, i.e., b and 9^. It is easy to show that the velocity at the 
edge of the plate reduces to = ±(ib/2) ( — I— e““‘®*). For the case of a 

120-degree camber this gives an absolute value of 0.433 with a velocity direction 
opposite to that expected at the edge of the plate ( Fig. 3.3). 
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Figure 3,3 Velocity profile along the radial line (nascent vortex on OM). 
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Placing the nascent vortex to the right of the radial line OM in the circle plane 
(or along the circular arc in the physical plane) always requires a stronger vortex to 
satisfy the Kutta condition and results in a tip velocity which is unrealistic both in 
magnitude and direction (Fig. 3.4). 

Placing the nascent vortex along the radial line OZ in the physical plane or 
outside the region .-X shifts the point of maximum velocity away from the edge of the 
plate (towards the downstream side). This, in turn, results in a leakage of fluid through 
the shear layer and requires a stronger vortex to satisfy the Kutta condition (Fig. 3.5). 

The entire region .A enclosed by the transformation of the radial line OZ in the 
physical plane and the radial line OM in the plane, is examined to determine the 
most appropriate positions of the nascent vortex. Figures 3.6 through 3.8 show a 
three-dimensional plot and the contour lines of F, L'(max), U(tip), and U(tip) as a 
function of the radial positions R and the angular positions R0 for those locations of 
the nascent vortex for which F< 1, U(max)'U(tip)<6. and U(tip) < 6. 

Figures 3.9 through 3.11 show three representative velocity profiles for three 
different positions of the nascent vortices in the region defined above. The most 
striking feature of these figures is that the maximum velocity near the tip can exceed 
considerably and unrealistically the velocity at the tip and that only for certain vortex 
positions does the maximum velocity (the velocity on the inner face of the camber at 
the tip) approaches smoothly the finite tip velocity. These calculations have shown 
that there is, in fact, a finite region in which the nascent vortices may be introduced in 
order to produce a tip velocity which is nearly equal to the maximum velocity in the 
vicinity of the tip. Clearly, it is only for unique combinations of the radial distance R 
and the angular position R0 that the said velocity ratio is equal to unity. 
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Figure 3.4 Velocity profile along the radial line (nascent vorte.x to the right of O.M). 



75 



-1000 -750 -500 -250 0 250 



nascont vortex 




GP^fl —3.15 



0.375 0.910 0.9-15 



0.380 1.015 

RE/RZ 



l.OSO 



1.085 




Figure 3.5 Velocity profile along the radial line (nascent vortex on OZ). 
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Figure 3.6 Three-dimensional plot of and the contour lines for F distribution. 
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Figure 3.7 Three-dimensional plot of the velocity ratio. 



7S 



.OJS l.U^U I.07S l.lUO 1.125 1.150 1.175 1.200 1.22S 1.250 





Figure 3.8 Three-dimensional plot of the tip velocity. 
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Figure 3.9 Velocity profile along the radial line (nascent vortex at c = 0.05, 60 = 25.0°). 
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.10 Velocity profile along the radial line (nascent vortex at C = 0.10, 50 = 25.0°). 
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Figure 3.11 Velocity profile along the radial line (nascent vortex at £ = 0.20, 56 = 25.0°). 



Figure 3.6 shows the required vortex strength as a function of a R and R0. Even 
though these calculations have been carried out with a single vortex, the subsequent 
calculations with larger number of vortices have shown that the relationship between 
the maximum velocity and the tip velocity remains practically invariant as long as the 
nascent vortices are always introduced at the fixed point which produces 
U(max) U(tip) ^ 1 for the single vortex. The reason for this is that the said velocity 
ratio is primarily dictated by the nascent vortex and the complex potential used in Eq. 
(3.29). to calculate the tip velocity, encompasses the effect of all the vortices in the flow 
field. In other words, it does not make anv difference whether the velocitv q .i.e.. (u 
+ iv^) in Eq. (3.20) is produced by the ambient velocity or by a large number of 
vortices in the field. 

The fact emerging from the foregoing analysis is that the nascent vortex cannot 
be placed arbitrarily (e.g., along the radial line in the circle plane or along the extension 
of the camber in the physical plane). Otherwise, the velocity distribution in the vicinity 
of the tip becomes unrealistic and unrepresentative of the evolution of the shear layers 
on either side of the camber. Furthermore, one is then forced to make arbitrary 
assumptions regarding the strength and the convection of the nascent vortices. 

The foregoing extensive analysis led to the conclusion that the nascent vortices 
should be introduced at r = 1.0925 and 0 = 0^ ± 2.08°, for the case of the 120-degree 
camber. To be more precise, the two nascent vortices are placed at the angular 
positions 0p = 0^^ — 2.08, and 0^ = 0^,-)-2.O8. during the period for which > 

0.2. For V, > Vp the positions of the two nascent vortices are switched to their 
corresponding images with respect to the radial line OM, i.e., they are placed at 0 . = 
0 . 2.08 and 0 . = 0 . - 2.08. The evolution of the verv earlv stages of the flow in 

51 qi s2 ^ 

the immediate vicinity of the tips of the camber is shown in Fig. 3.12. 

In Fig. 3.12a the velocity field is a consequence of the first two nascent vortices 
introduced at the points noted above. Figures 3.12b through 3.1 2d show the 
development of the flow field and the starting tip vortex subsequent to the introduction 
of the 4^, 8*, and 11* nascent vortex, respectively. 

There is not a unique procedure for relating the rate at which vorticity is shed 
into the wake, the Kutta condition, the velocitv with which the nascent vortices are 
convected, and the time interval for the convection, all of which help to simulate the 
e.xperimentally observed features of the free shear layers. Fage and Johansen (1928). 
through quite ingenious experiments with steady flow about various bluff bodies, have 



83 




Figure 3.12 Velocity field in the vicinity of the tip. 
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^0.48 T=064 



shown that vorticity is shed from the two sides of an axisymmetric body (a circular 
cylinder) or a sharp-edged body (a plate normal to the flow) at the same rate; that the 
motion in a sheet is steady near the body, except possibly near the inner edge of the 
shear layers; that fluid flows into a sheet through both edges, but at a greater rate 
through the outer edge; that at each section of the sheet the velocity rises from a small 
value to a well-marked maximum value (approximately Vj U = 1.45) and then very 
slowly decreases to about 1.35 within a distance of approximately y~2c, where the 
breadth of the sheet reaches a value of A ~ c; and finally, that the velocity Vj at the 
outer edge of the sheet is much larger than the velocity at the inner edge (except 
during the deceleration period of the flow) and V," may be ignored in calculating the 
vorticity flux from cF ct =0.5 ( Vj- — \%-). 

In the present calculation the vortex strength, the velocities on either side of the 
shear layer, and the time interval are related by 

r = 0.5 (V|- - Vr)-At (3.37) 



in which F is the strength of the nascent vortex, Vj = U(tip) and V, is the velocity at 
the downstream face of the camber near the tip. The velocity V, can be calculated 
correctly in a number of ways, to be described later. Suffice it to note that in general 
\\ is very small (for steady flow) and that the method of its calculation has ver>' little 
or no influence on the strength of the nascent vortex or on the time interval to be used 
for a given vortex strength. The velocity V, becomes important only when the wake 
begins to move towards the camber (i.e.. during the period of flow deceleration). 

To explain the computational details of the method let us consider a particular 
time t after the start of the motion and assume t to be sufficiently large so that there 
are a number of vortices in the wake. Then the appearance and convection of the 
vortices proceed as follows: 

(1) Determine the strength of the nascent vortices from the Kutta condition ,Eq. 
(3.20), in which Cq is a known fixed position for each nascent vortex ; 

(2) Place the nascent vortices at Cgp and Cg^ and calculate the velocity Vj at the 
two edges of the plate ; 

(3) Calculate V.,, representing the velocity at the inner boundary of the shear 
layer, as the average of the velocities at three points along the radial line OZ 
in the physical plane i.e., at r = 1, 1.05, and 1.1; 
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(4) Calculate the time interval, for each edge of the plate, from Eq. (3.3"! using 
the known values of F and the velocities and V.,. Store the average of the 
two time intervals for use in subsequent calculations: 

(5) Calculate the velocity induced at the center of all other vortices; 

(6) Convect the two nascent vortices with a velocity 0.5(V^ ^ V., ) for an average 
time interval At (note that the vorticity is convected with the average velocity 
of the shear layer). If the distance travelled by a nascent vorte.x is not within 
0.05 = 0.01. it is convected twice for a time interval At 2. The subsequent 
convections of the nascent vortices are miade using the velocity induced at 
their center. 

(7) Convect all other vortices for the same time-interval At using a second order 
schem.e given by 

z(t-^At) = z(t) 0.5 [3z(t) - z(t-At)]«At (3.38) 

in which z = ti iv. 

(S) Remove the vortices fro.m the calculation whenever they come nearer than 
0.05 to the camber in the physical plane (except the first 20 vortices from the 
tip); 

(9) Coalesce the same sign vortices with a separation of less than 0.05 (in the 
physical plane, except the first 20 vortices); 

(10) Calculate the tangential velocities and pressures on the inner and outer faces 
of the camber. Determine the drag and lift forces through the integration of 
pressure and through the use of the rate of change of impulse. .Make plots of 
suitable variables (e.g.. velocity distribution near the tip. variation of nascent 
vortex circulation with time, evolution of the wake, etc.V. 

(11) Check the flow conditions to determine the state of the calculations: 

(aj If Vj - V, > 0.2 repeat the foregoing steps; 

(b) Stop the introduction of nascent vortices if 0 < - V, < 0.2 and 

return to step No. 5; 

(c) If V, > switch the angular positions of the nascent vortices to their 
image points. Calculate V, as the average of the three velocities, at the 
upstream, side of the tip of the camber, at three radial locations (0.95. 
0.9. and 0.S5) and repeat the foregoing steps; and 



S6 



(12) Make plots of the variations of various flow parameters (e.g., tip velocity, 
nascent vortex circulation, evolution of the wake, force coefficients, etc.) and 
terminate the run. 

The foregoing steps are quite general and can be used for any camber, provided 
that the optimum points of placement of the nascent vortices are determined through a 
similar analysis for the desired camber angle. 
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IV. DISCUSSION OF RESULTS 



A. NUMERICAL AND PHYSICAL EXPERIMENTS 

The calculations were carried out for a time-dependent normalized velocity given 
by 




for 



T"'= < 8.65 

c 



and 



(4.1) 




0. 1539(T’'’ - 8.65) + 0.0053 1(T=' - 8.65)‘ 



(4.2) 



which corresponds to that encountered in a series of experiments carried out in a 
vertical water tunnel. A detailed description of the equipment and procedures is given 
by Sarpkaya and Ihrig (1986) and will not be repeated here. Evidently, the calculations 
can be carried out for any specified variation of the velocity. For the case under 
consideration, the flow begins to decelerate at T* = 8.65 and the velocity of the 
ambient flow reduces to zero at about T* = 19, (see Fig. 4.1) 

The com.puter program provided, at times specified, the positions of all the 
vortices, the rate of shedding of vorticity from the tips of the camber, the velocity 
distribution on the upstream and downstream faces of the camber, the total and 
differential pressure distributions, and the force coefficients. 

Figures 4.2 through 4.4 show, at T* = L’^t 'c = 4.35, the evolution of the wake, 
the tangential velocities at the upstream and downstream faces of the camber, the 
velocity profile along the radial line passing through the tip. and the total and 
differential pressure distributions (these plots are available at every time step but are 
not reproduced here for sake of brevity ). These and other figures show that the 
characteristics of the flow develop symm.etrically prior to the onset of deceleration (T* 
< 8.65) and the diflerential pressure is positive evervxvhere (i.e.. the pressure inside the 
camber is larger than that outside). 
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Figure 4.1 Variations of the velocity and the acceleration. 
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Figure 4.2 Position of vortices and tangential velocity distribution. 
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Figure 4.3 Velocity profile along the radial line passing through the tip. 
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Figure 4.4 Pressure distribution on the upstream and downstream faces. 
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Following the onset of deceleration (see e.g., Figs. 4.5 through 4.7 at T* = 
10.S4\ the difTerential pressure near the axis of the camber becomes increasingly 
negative. The reason for this is that the deceleration of the flow brings the vortices 
closer to the camber. The significance of this result is that had the model been flexible 
(as in the case of a parachute) the central part of the camber would have collapsed as a 
result of the particular deceleration it is subjected to. Evidently the collapse 
phenomenon would not have remained symmetrical, as evidenced by field experiments 
with axisymmetric flexible parachute. Furthermore, it would have required the analysis 
of the flow about a flexible camber. It is because of this reason that in the present 
analysis the camber is assumed to remain rigid. 

For T* larger than about 13 (for the ambient flow under consideration), the 
velocities induced at the downstream edges of the camber by the large vortices moving 
sideways and towards the camber give rise to oppositely-signed vorticity. This, in turn, 
leads to the rapid growth of the secondary vortices (see e.g., Figs. 4.S through 4.10 at 
T* = 17.86). The secondar\' vortices are relatively weaker than the priman.' vortices 
partly because they have been in existence only for a short time and partly because the 
vorticity flux is not as large as that in the primary shear layers. Consequently , the 
centroid of the secondary' vortices tends to orbit about the centroid of the primary 
vortices. 

A comparison of Figs. 4.7 and 4.10 shows that the region of negative differential 
pressure grows with time and occupies a large central portion of the camber. In fact, 
the drag force acting on the camber becomes negative, as it will be seen shortly. 

Figures 4.11 through 4.16 show, at suitable times, the velocity field about the 
camber. The rapid growth of the wake during the period of steady uniform flow is 
exhibited in Figs. 4.11 and 4.12. Figure 4.13 nearly corresponds to the time at which 
the deceleration is imposed on the flow. Figures 4.14 through 4.16 show clearly the 
backward motion of the primary’ vortices and the rapid growth of the secondar.- 
vortices. It is seen from Fig. 4.16 that the fluid motion is entirely due to the motion of 
the vortices in the flow" field. The two vortices on each side of the camber form a 
counter-rotating couple and remove themselves rapidly from the field under the 
influence of their mutual induction velocity. Subsequently, the absolute value of the 
the difTerential pressure begins to decrease. Eventually, the differential pressure reduces 
to zero even.'where on the camber as the conditions approach to that of a body in a 
fluid at rest. 
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Figure 4.5 Tangential velocity at T* = 10.84. 
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Figure 4.6 Tip velocity at T* = 10.84. 
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Figure 4.7 Pressure distribution at T* = 10.84, 
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Figure 4.8 Tangential velocity at T* = 17.86. 
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Figure 4.9 Tip velocity at T* = 17.86. 
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Figure 4.10 Pressure distribution at T* = 17.86. 
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Figure 4.11 The velocity field about the camber at T* = 4.35. 
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Figure 4.12 The velocity field about the camber at = 6.05. 
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The remainder of the discussion of the results will deal with the variation with 
time of the tip velocities and the drag and lift forces. 

Figures 4.17 and 4.18 show the velocities and V, as a function of T’''. The tip 
velocity V'j decreases from an initially large value of about 3.5 to a nearly constant 
value of about 1.5 just prior to the onset of deceleration. Subsequently, decreases 
rapidly during the period of deceleration and prior to the inception of the secondarv' 
separation. Then increases to about 2 because of the backward motion of the large 
vortices near the tips of the camber. Finally, Vj decreases once again as the primary 
and secondary vortices move sideways and away from the tips of the camber due to 
their mutual induction (see Fig. 4.16). 

The variation of with T* is significant only during two, relatively short, time 
intervals: at the start of the motion and at the start of the deceleration. These are the 
periods during which the vorticity flux changes rapidly in order to maintain the Kutta 
condition. During the remainder of time V, is negligibly small, as expected on the 
basis of the pioneering e.xperiments of Fage and Johansen (1928) with steady flow over 
various types of bluff bodies. 

Figures 4.19 and 4.20 show the variation of the drag and lift coefficients as a 
function of time. The former is based on the integration of pressure and the latter on 
the rate of change of impulse. The drag coefiTicient calculated through the use of the 
rate of change of impulse is somewhat larger than that obtained through the 
integration of the instantaneous dilTerential-pressure distribution. This is due to the 
fact that the impulse expression includes the rate of change of circulation between two 
successive time steps whereas the pressure expression does not. It is a well-known fact 
that in real fluids the memory of the fluid resides in its vorticity. Whereas in inviscid 
fiows there is no memory and the dynamic characteristics of the How (pressures and 
forces) are functions of only the instantaneous state of the flow. The analysis 
presented herein is for an inviscid fluid even though the phenomenon concerns the 
motion of a real fluid. The question of whether the rate of change of circulation 
should be included or excluded in the discrete vortex analysis (first discussed by 
Sarpkaya in 1968) is an unsettled issue. It appears that only the comparisons with 
experiments can clarify the question. 

Figures 4.19 and 4.20 also show that rises rapidly (due to the rapid 
accumulation of vorticity in the growing vortices) and begins to decrease as the 
vortices develop under the inlluence of a constant ambient velocity. Then the force 
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Figure 4.13 The velocity field about the camber at T* = 8.55. 
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Figure 4. 14 The velocity field about the camber at T* = 14.20. 
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Figure 4.15 The velocity field about the camber at T* 
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Figure 4.16 The velocity field about the camber at T’^ = 21.50. 
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Figure 4.17 Variation of with time. 
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Figure 4.18 Variation of with time. 
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decreases sharply at the onset of deceleration and goes through zero near the middle of 
the deceleration period (T* = 11). The force acquires its largest negative value 
towards the end of the deceleration period. Subsequently, the force gradually decreases 
to zero. 

Also shown in Figs. 4.19 and 4.20 is the variation of the lift force. It is negligible 
even in the later stages of the motion. This is primarily due to the fact that there is 
not sulTicient time for the development of alternate vortex shedding either during the 
period of steady flow or during the period of rapid deceleration. 

Figure 4.21 shows a comparison of the calculated (through pressure integration) 
and measured drag coefficients. In general the agreement between the calculated and 
the measured drag coefficient is quite good. In the time intervals between 13 and 16 
and between 19 and 22. the calculated is somewhat larger. The reason for this is as 
follows. In the said time intervals the drag coefficient is relatively small and the 
viscous effects are relatively important in dissipating the vortices. This is not taken 
into consideration in the numerical analysis. It is possible to introduce a small 
artificial reduction in circulation in order to bring the calculated and measured values 
into closer agreement. This has been avoided in the present analysis in order to keep 
the discrete vortex analysis as pure and simple as possible. Figure 4.21 also shows 
that the calculation of the drag coefficient through the integration of pressure is 
superior to that through the use of the rate of change of impulse. 

Finally, a comparison is made between the calculated and photographed flow 
fields at corresponding times. Figures 4.22 through 4.24 show at times T* = 6.05, S.55 
and 16.30 the flow in the immediate vicinity of the camber (plotted to the same scale), 
shows that the agreement between the calculated and observed flow fields is indeed 
ver>’ good. 

B. CONCLUDING REMARKS 

The results presented in this section have shown that the discrete vortex model 
can be used with confidence to predict the evolution of the wake about a cambered 
plate immersed in a an arbitrary time-dependent flow. The evolution of the wake is 
remarkably similar (including the formation of the secondary vortices) to that obtained 
in flow visualization experiments. The drag coefficients resulting from the analysis and 
experiments agree reasonably well. This agreement can be improved with the 
introduction of a small circulation dissipation. The drag coefficient calculated through 
the integration of the instantaneous pressure distribution agrees more closely with that 



IU9 




Figure 4.19 Drag coefficient calculated from pressure integration. 
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Figure 4.20 Drag coefficient calculated from impulse expression. 
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Figure 4.21 Calculated and measured drag coefTicients. 
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Figure 4.22 Comparison of flow fields at T* - 6.05. 
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Figure 4.23 Comparison of flow fields at T’* — 8.55. 
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Figure 4.24 Comparison of flow fields at T* - 16.30. 
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obtained experimentally and points out the fact that the inclusion of the rate of change 
of circulation term in the impulse method is not in conformity with the behavior of the 
How. 

The development of negative diflerential pressures near the central region of the 
camber is thought to be primarily responsible for the inception of the partial collapse 
of a parachute at high rates of deceleration. This phenomenon takes place even when 
the total drag force acting on the parachute is still positive. The sample analysis 
presented herein also shows that the negative differential pressure can cover a large 
region of the parachute and even result in negative drag. The basic idea emerging from 
the analysis reported herein is that the designs which incorporate into them the idea of 
delaying or preventing the return of the shed vortices to the canopy (e.g. .porosity 
management, change of deceleration history, parachute shape, dissipation and or 
destruction of the organized wake) will be the ones which could avoid the collapse 
phenomenon. Extensive analysis and small scale experiments coupled with few 
judiciously selected field tests may help to arrive at practically and phenomenologically 
sound parachute designs. 
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V. CONCLUSIONS 



In this chapter the conclusions reached in the preceding sections are summarized 
and the areas in which the greatest need for further work is perceived are pointed out. 

.‘\ brief review has shown that two-dimensional finite difference and finite element 
methods can simulate successfully the low Reynolds number flows but serious 
impediments remain to extending them to high Reynolds number range that is of 
practical interest. The major obstacles are: the large computational requirements 
which grow rapidly with Reynolds number, stability problems, difi'iculty of 
implementing boundary conditions and artificial viscosity. 

The discrete vortex model is seen to avoid many of the problems of the Eulerian 
finite difference methods but are subject to some of their own. Chief among these is 
the introduction of vorticity, determination of the separation points, excessive 
computation time for the convection of vortices, and the need to introduce artificial 
dissipation to bring the measured and calculated results into closer agreement. 

The first section of this study dealt with the numerical simulation of a 
sinusoidally oscillating flow about a circular cylinder. 

The discrete vortex model has been combined with the boundary layer 
calculations and the positions of the separation and stagnation points have been 
calculated as accurately as possible for a Keulegan-Carpenter number of K = 10. The 
results have accurately predicted the formation of a half Karman vortex street in the 
transverse direction. The calculated positions of the vortices were found to be in good 
agreement with those obtained experimentally. The measured and calculated drag 
force and the differential pressure distributions showed reasonably good agreement. 
The results have also shown that the effect of the backward convection of a large 
vortex over one side of the cylinder is indeed very pronounced on all the measured and 
computed characteristics of the flow. This is one of the most important reasons as to 
why the .Morison's equation (see e.g.. Sarpkaya & Isaacson 1981) fails to represent the 
in-line force acting on the cylinder with reasonable accuracy particularly in the range S 
< K < 13. 

The number of numerical parameters involved in the simulation made a complete 
parametric study of their effects impossible. However, several tests were made in which 
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the time step and convection scheme were varied. The results showed e.xpected 
improvements as the time step was reduced and the order of convection scheme was 
increased but no undue sensitivity was observed. 

The application of the discrete vqrte.x model to a particular time-dependent flow 
past a camber required an extensive 'study of the velocity field in the vicinity of the 
sharp edges of the camber. The results have shown that the nascent vortices can be 
introduced only at judiciously selected points. The numerical experiment predicted 
satisfactorily the evolution of the wake and the forces acting on the body. 
Furthermore, the calculations have provided a plausible explanation for the cause of 
parachute collapse, a phenomenon which has provided the impetus for the simulation 
described herein. 

In the course of the present work, two aspects of the discrete vortex model were 
identified that seem to be in greatest need for further investigation. The first is the 
determination of the mobile separation points on a body without sharp edges. The 
existing methods cannot be considered satisfactoiy for the prediction of separation in 
unsteady laminar and or turbulent boundary layers. The use of Pohlhausen's method 
seemed to produce good results (as judged by the experimental data), but it is the most 
difficult part of the model to defend. 

The second aspect of the model which requires further work is the discovery of a 
systematic and conceptually satisfactory method to reduce vorticity. It seems clear, at 
least on the basis of experiments, that there is a physical basis for dissipating vorticity. 
but it remains to be proved that the global error involved in doing so is reasonable and 
in conformity with the behavior of nature. Future efforts will probably have to be 
preceded by careful measurements of the finer details of the unsteady flow field. 
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APPENDIX A 

NUMERICAL RESULTS 



This section presents additional figures showing the evolution of sinusoidally 
oscillating flow about a cylinder. It is seen that the convection of the previously shed 
vorte.x over the shoulder of the cylinder precipitates earlier separation, establishes a 
preferred position for the ne.\t dominant vorte.x, and gives rise to additional primary 
and secondarx’ separation points. 
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Figure A.l Position of vortices, velocity and pressure distribution at T* 1.025. 
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Figure A. 2 Position of vortices, velocity and pressure distribution at T*= 1.100. 
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Figure A. 3 Position of vortices, velocity and pressure distribution at T* l.2oO. 
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Figure A. 4 Position of vortices, velocity and pressure distribution at T*= 1.350. 
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Figure A. 5 Position of vortices, velocity and pressure distribution ar T*- 1.525. 



124 



PTOTfil. VEL Y/C 




' -12 -io -3 -5 -2 0 2 4 S 

x/c 



1.6QGO 



I ^ n 

3 10 



1 

12 



TRNGENTIRL VELOCITY 




TOTRL PRESSURE 




Figure A. 6 Position of vortices, velocity and pressure distribution at T*= 1.600. 
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Figure A. 7 DiiTerential pressure distributions at T*= 1.350 and 1.600. 
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Figure A. 8 Position of vortices, velocity and pressure distribution at T*- 1.756. 
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Figure A. 9 Position of vortices, velocity and pressure distribution at T’” - l.SoO. 
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Figure A. 10 Position of vortices, velocity and pressure distribution at T*= 1.950. 
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Figure A. 11 Position of vortices, velocity and pressure distribution at T* = 2.250. 
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Figure A. 12 Position of vortices, velocity and pressure distribution at T* = 2.300. 
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Figure A. 13 Position of vortices, velocity and pressure distribution at T* = 2.450. 
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Figure A. 14 Position of vortices, velocity and pressure distribution at T* = 2.550. 
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Figure A. 15 Position of vortices, velocity and pressure distribution at T*- 2.725. 
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Figure A. 16 Position of vortices, velocity and pressure distribution at T* = 2.925. 
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Figure A. 17 Position of vortices, velocity and pressure distribution at T’' = 2.975. 
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